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ABSTRAGT 

Photometry alone is not sufficient to unambiguously distinguish between ultra-faint star clusters 
and dwarf galaxies because of their overlap in morphological properties. Accurate measurements of 
the intrinsic velocity and metallicity dispersions are required to robustly classify such objects. Here 
we report on VLT/GIRAEEE spectra of candidate member stars in two recently discovered ultra¬ 
faint satellites Reticulum 2 and Horologium I, obtained as part of the ongoing Gaia-ESO Survey. We 
identify 18 members in Reticulum 2 and 5 in Horologium I. We find Reticulum 2 to have a velocity 
dispersion of 3.22 /q 49 kms“^, implying a M/L ratio of ^500. We have inferred stellar parameters 
(Teff, log^, [Ee/H], [g/Ec]) for all candidates and we find Reticulum 2 to have a mean metallicity of 
[Fe/H] = -2.46tno, with an intrinsic dispersion of 0 . 29/9 ^5 dex, and is G-enhanced to the level of 
[g/Ec] ^ 0.4 dex. We conclude that Reticulum 2 is a dwarf galaxy. We also report on the serendipitous 
discovery of four stars in a previously unknown stellar substructure near Reticulum 2 with [Ee/H] ^ — 2 
and Vhei rsj 220 km s far from the systemic velocity of Reticulum 2 at Vhei = 64.7/J;gkm s Eor 

Horologium I we infer a velocity dispersion of a (U) = d.Q/gig kms~^ and a consequent M/L ratio of 
^600, leading us to conclude that Horologium I is also a dwarf galaxy. Horologium I is slightly more 
metal-poor than Reticulum 2 ([Ee/H] = —2.76/o.i) arid is similarly G-enhanced: [g/Ec] ^ 0.3 dex. De¬ 
spite a large error-bar, we also measure a significant spread of metallicities (cr ([Ee/H]) = O.I 7 /o.o 3 
which strengthen the evidence that Horologium I is indeed a dwarf galaxy. The line-of-sight velocity 
of Reticulum 2 is offset by some 100 kms“^ from the prediction of the orbital velocity of the LMC, 
thus making its association with the Cloud uncertain. However, at the location of Horologium 1 , 
both the backward integrated LMC’s orbit and the LMC’s halo are predicted to have radial velocities 
similar to that of the dwarf. Therefore, it is very likely that Horologium 1 is or once was a member 
of the Magellanic Eamily. 

Subject headings: Galaxy: halo, galaxies: dwarf, globular clusters: general, galaxies: kinematics and 
dynamics, stars: abundances, galaxies: abundances 


1. INTRODUCTION 

The total inventory of satellites associated with the 
Milky Way remains incomplete. This is particularly 
true for the faintest systems, as observations are inher¬ 
ently biased towards finding and characterising intrinsi¬ 
cally bright satellites. As a consequence, any inferences 
on the total satellite luminosity or mass distributions 
strongly depend on systematic and selection effects in 
the least luminous systems. Deep, uniform photometry 
is required to find these relics, and spectroscopy is re¬ 
quired for proper characterisation. 


Wide-field photometric surveys can b e extremely suc- 
cessful at finding Galactic sat ellites dWillman ^t ah 
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data has been used by two i ndependent g:rQiips to find 


at least another 9 satellites 1 

jKoposov et al. [2015 

The 

DES Collaboration et al.[[201. 

5[). T'he wide-field imaging 


number of smaller-scale surveys of the Milky Way halo, 
increasing the tally of Galactic satellites at a breakneck 


al. 

2Ul5l. 

abl 

y feet). 


Martin et 


Many of these latest discoveries are remark- 


earning the name of ‘ultra-faint’ satellites. 

As more ultra-faint satellites (UFS) have being discov¬ 
ered, complications have arisen in trying to accurately 
classify them. A tenuous overlap between the effective 
radii and absolute magnitudes of faint globular clusters 
and ultra-faint dwarf galaxies has emerged. It is now 
crystal clear, the morphology of ultra-faint systems near 
the valley of a mbiguity cannot be c lassified from pho¬ 


tometry alone (Gilmore et al. 2007). Kinematics and 


chemistry are required to distinguish between these clas 
sifications. A large dispersion in overall metallicity is rep 
resentative of extended star formation in a dwarf galaxy- 
like environment that is massive enough to retain super¬ 
nova ejecta, providing a key diagnostic for di stinguish¬ 
ing globular clu sters and dwarf galaxies (e.g., Willman 


& Strader 2012). Indeed, spectroscopy is essential for a 
large number of confirmed members in order to precisely 
measure velocity and chemical dispersions, estimate the 
dark matter content, and explore the star formation his¬ 


tories of these ancie nt systems (e.g., [Kirby et aL 
Tollerud et al.|[2M^ ) 


2011 


The UI^'Ss have been the focus of attention of Galactic 
archaeologists world wide for less than a decade. Wor- 
ryingly, during this short history, many of their spec¬ 
troscopically determined properties have continued to 
evolve. For example, the early studies of the Bootes 1 
dwarf spheroid al reported a yeloci t v dispersion as high 


as 6.5 km s ^ (^Muhoz et al. 


2006 


Ho wever, an independent and novel study by [Koposov et 


Martin et al. 2007). 


al. (2011) revealed that the dwarf’s internal kinematics 


IS potentially dominated by a stellar population with a 
velocity dispersion as low as 2.4 km s~^. Similarly, for 


the Segue 2 satellite, Belokurov et al. (2009) gave an estl 
mate of 3.4 kms“^ and a warning of a potential contam- 
ination from the s urrounding Tri-And stream. In fact, 
Kirby et al. (2013) later showed the velocity of Segue 2 
IS consistent with zero, thus ruling out the presence of 
any significant amount of dark matter in the system. 
Undoubtedly, robust uncertainties on individual stellar 
velocity measurements are paramount to accurate char¬ 
acterization of the kinematics of these systems. Further¬ 
more, inferences are susceptible to low-number statistics, 
as well as contamination from foreground stars or binary 
systems. 


The Gaia -ESO Survey (Gilmore et al. 2012 Randich 


et al. 2013) has been obtaining high-resolution spectra 


using the Very Large Telescope (VLT) in Chile since Jan¬ 
uary 2012. The primary scientific goal of the Gaia-ESO 
Survey is to characterise the major constituents of the 
Galaxy, and to understand these components in the con¬ 
text of the Milky Way’s formation history. To that end, 
more than 10^ Milky Way stars are homogeneously tar¬ 
geted as part of the Gaia-ESO Survey, including all ma¬ 
jor structural components: open and globular clusters. 


the disk, bulge, and the halo. Understanding the role of 
ultra-faint satellites is indeed important in this context, 
as they inform us of star formation in isolated environ¬ 
ments as well as the accretion history of the Galaxy. 

Here we report on Gaia-ESO Survey observations of 
two of the recently discovered ultra-faint dwarf satel¬ 
lites: Reticulum 2 and Horologium 1. Reticulum 2 is a 
mere ^30 kpc away towards the Large Magellanic Cloud 
(LMC) with a My = —2.7 ± 0.1, and is suspected to be 
very dark matter-dominated. Indeed, Reticulum 2 is of 
particular scientific interest given the recently reported 
detection o f dark matter annihilation (jGeringer-Sameth 


et al. 


2015 The Eermi-LAT Collaboration et aT 2015” 
Horologium 1, also located towards the LIVIC, is more dis- 
tant at 79 kpc but given somewhat brighter luminosity of 
My = —3.4±0.1 and a visible giant branch, spectroscopy 
of candidate red giant branch (RGB) stars is accessible 
from the VLT in a standard Gaia-ESO Survey observing 
block. While our primary goal is to establish the true na¬ 
ture of these faint stellar systems by gauging the amount 
dark matter they might contain, we also intend to test 
the hypothesis that the dwarfs have once been part of 
the Magellanic group. 

We outline the target selection and the data that we 
will subsequently analyse in Section A detailed de¬ 
scription of our analysis in outlined SectionWe discuss 
the interpretations of our results in Section]^ 


2. OBSERVATIONS 


The data were obta ined in 0.8^^ seeing using the 
EL AMES spectrograph (jPasquini et al.||2002|) on the 8.2 
m Very Large Telescope (VLT'/UT'2) at Gerro Paranal, 
Chile. Candidate members of both satellites were tar- 
getted using otherwise unallocated Milky Way fibres on 
Eebruary 6 and March 8-10, 2015 as part of the on- 
going Gaia-ESO Su r vey (ESQ Propamm e 188.B-3002 


Gilmore et al. 


2012 


2013). In the field 


^ [Randich et al 
GESJy[W_033542_540254 109 fibers were allocated to sci 
ence targets, with 25 allocated to Reticulum 2 can¬ 
didates, and the rest allocated to standard Gaia-ESO 
Milky Way targets. In the field GES_MW_025532_540711, 
the total number of allocated fibers was 107, with 
18 fibers assigned to Horologium 1 candidates. The 
HRIO and HR21 setups were ermloyed, providing high- 
resolution {IZ ~ 19800 and 1620Cn) spectra in wavelength 
regions of 5334-5611 A and 8475-8982 A, respectively. 

The candidate satellite members were selected using 
a broad color-magnitude mask based o n the best fit¬ 
ting is ochrone and distance modulus from [Koposov et ah 
(2015[) (hereafter K15). We also required that the targets 
were located within 10 - 15' on the sky from the satellite 
center. Eigure[^ shows the color-magnitude distribution 
of stars near the center of both systems. Candidates that 
were observed spectroscopically are highlighted, as well 
as those which we later confirmed to be members. 

The data were reduced using standard procedures per¬ 
formed for all other Gaia-ESO Survey GIRAEEE ob¬ 
servations. This process includes bias correction, flat- 
fielding, object extraction, sky subtraction, scattered 

^ Note however that the resolving power 
and sensitivity of the GIRAFFE instrument has 

been recently improved due to refocusing, see 

http://eso.org/sci/pubhcations/announcements/sciannI50I3.html 
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Reticulum 2 Horologium 1 



Fig. 1. — Color-magnitude distribution of stars near the centers 
of Reticulum 2 and Horologium 1 satellites. Large black circles 
indicate candidates that were observed with VLT/GIRAFFE, and 
the red symbols are those we later confirmed to be member s of each 
system. Eac h panel shows a 12.5 Gyr PARSEG isochrone |Bressan| 
|et al.| (|2012|) with a [Fe/H] = — 2 that has been shifted according 
to the distance modulus from K15. 


light correction, and wavelength calibration. The spec¬ 
tra are then corrected for barycentric mo tion and re¬ 
samp led onto a common wavelength scale (Lewis et ah 
pMst - Generally the mean signal-to-noise (b/M) ratio 
per pixel is quite low: 5/13 for Reticulum 2 candidates in 
HR10/HR21 respectively, and 4/10 for the Horologium 1 
candidates. However the brig htest confirmed members 
of either satellite (Section |3.3[ ) have S/N ratios of 26/62 
(Reti 4jHR10/HR21) and 10/28 per pixel (Horo 10; see 
Figure 1^. 

3. ANALYSIS 
3.1. Kinematics 

The velocity analysis has been performed using stan¬ 
dard Gaia-ESO Survey radial velocity (RV) pipeline. Al- 
th ongh the ful l detai ls of this pipeline will be presented 
in Lewis et al. (|2015 ), here we outline the principle com- 


ponents. 'I'he algorithm is based on a direct pix el-fitting 
procedure implemented by Koposov et al. (2011) (see also 


Koleva et al. 2009 Walker et al.'2015a) , and e mploys 
the FHUE NIX library of model stellar spectra (jHusser 
et al. All of the observed spectra are fitted (using 

a metric) by model templates that are interpolated 
from the spectral library, and the continuum is modelled 
by a high-degree polynomial. The maximum likelihood 
point is foun d by using Nelder-Mead algorithm (Nelder & 
Mead |1965), using several starting points to avoid being 
trapped in local maxima. The location of the maximum 
likelihood (y^ minimum) and the Hessian of the likeli¬ 
hood surface are used to provide the best-fit velocity, 
errors and estimates of the stellar atmospheric parame¬ 
ters. 

The crucial ingredient in correctly extracting the kine¬ 
matics of ultra-faint satellites is a proper understa nding 
of the uncertainties of RV meas urements (see e.g., Geha 
et al.|2QQ9 ; Koposov et al.|2QlT ). Because the Gaia-EbO 


Survey has already observed many thousands of stars, 
some of them repeatedly and with different instrument 
configurations, our experience with the Survey provide 
us with a very good understanding of the RV precision 
achievable. Eor this paper, however, we focus only on 
data in just two Gaia-ESO Milky Way fields: those with 
Horologium 1 and Reticulum 2 candidates. 

In addition to the standard data processing steps per¬ 
formed for the Gaia-ESO Survey, there are three subtle 
points that are important for this study: 

• Spectral covariance. The standard Gaia-ESO re¬ 
duction pipeline rebins the spectra to a common 
wavelength mapping with a fixed step size. Eor 
HRIO the common wavelength scale extends from 
5334-5611A and for HR21 the boundaries are 8475- 
8982A, both with a spacing of 0.05A. While con¬ 
venient for some analyses, the rebinning procedure 
introduces correlated noise/covariance in the spec¬ 
tra and reduces the effective information content of 
the spectra. Eor example, the rebinned HR21 spec¬ 
tra has 10141 pixels, the rebinned HRIO spectra 
has 5541 pixels, while the original spectra are just 
4096 pixels. We can account for this correlation 
by modelling the spectra with the full covariance 
matrix, or approximate it. Our tests found that if 
we fit the spectra using the full covariance matrix 
of the data and the posterior/likelihood is properly 
behaved (e.g., uni-modal and close to a Gaussian), 
then the effect of pixel covariance is equivalent to 
scaling the errors by a fixed constant: 1.5 for HRIO, 
2.0 for HR21. These numbers are approximately 
equal to the ratio of rebinned and original pixels. 
We adopt this scaling throughout the rest of our 
analysis. See the end of this section for the verifi¬ 
cation of the results. 

• Systematic error floor. It is well-known that al¬ 
though the formal RV precision derived from cross¬ 
correlation or pixel-fitting methods can be almost 
arbitrarily small for sufficiently high S/N spectra, 
the actual precision achievable with most spectro¬ 
graphs is generally limited by systematic effects. 
This includes instrument flextures, uncertainties 
in the wavelength calibrations. Line Spread Eunc- 
tion (LSE) variation/asymmetry and template mis¬ 
matches. This systematic component has to be in¬ 
cluded in the total error budget. We have found 
this systematic error to be around 300 m s“^ from 
large numbers of Gaia-ESO Milky Way spectra. It 
is important to note that this systematic compo¬ 
nent is not expected to be present when comparing 
RVs obtained from spectra using the same setup 
in sequential exposures, but it becomes important 
when comparing RVs from different nights, or be¬ 
tween HRIO and HR21 setups. We include this 
systematic error floor in suitable comparisons here¬ 
after. 

• RV offset between HRIO and HR21. Over the 
course of the Gaia-ESO Survey, it has been dis¬ 
covered that there is a small systematic offset of 
400 m s“^ between the RVs measured in the HRIO 
and HR21 setups. The cause of this offset is not 
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Fig. 2. — An example showing the observed (black) and maximum a posteriori model (red) spectra for a confirmed Reticulum 2 (top) 
and Horologium 1 (bottom) member. These spectra have the highest S/N ratio of confirmed members in each system. 


well-established yet. This correction was applied 
to the radial velocities (the HR21 velocities have 
been shifted by —400ni s“^). 

After applying the aforementioned corrections, we can 
confirm whether the RVs measured in repeated expo¬ 
sures match within the precision quoted by our error 
bars. To test this we have collated all the spectra in 
the GES_MW_033542_540254 and GES_MW_025532_540711 
fields (i.e., including both standard Gaia-ESO targets 
and possible satellite member stars). The top panel of 
Figure shows the distribution of velocity differences 
{V 1 — V 2 ) scaled by the RV error (^cr(Vi)^ + cr(V 2 )^) for 
repeated HRIO exposures. The middle panel of the figure 
shows the same for the HR21 setup. The bottom panel 
shows the distribution of normalised velocity differences 
for HRIO exposures versus HR21. In all the panels the 
red curve shows a standard normal distribution with zero 
mean and unit variance. In all cases the distributions are 
indeed well described by Gaussians, confirming that our 
error model provides a correct description of the velocity 
uncertainties. 

The final radial velocities for all the observed Reticu¬ 
lum 2 and Horologium 1 candidate members are provided 
in Table and refer to the weighted means of the HRIO 
and HR21 measurements and take into account all the 
error-terms mentioned above. 


3.2. Chemistry 

We used a generative model to infer the stellar 
parameters for all stars. The model is described 
as follows. For a given set of stellar parameters 
uj = {Teff,log^, [Fe/H], [of/Fe]} we first produce a fiux- 
normalised synthetic spectrum at wavelengths A 

by interpolating spectra from a surrounding grid. The 
synthetic spectra were calculat ed as per the AMBRE 
grid (see de Laverny et al.||2012 for details). This high- 
resolution (7^ > 300, 000) grid was synthesised spe cifi- 


cally for the Gaia-ESO Survey using Turbospectrum (Al- 



Fig. 3.— The results of radial velocity precision tests done using 
repeated observations. The three panels show the distribution of 
the RV difference normalised by their errors as measured in the 
HRIO (top panel), HR21 (middle panel) and HR10-HR21 config¬ 
urations (bottom panel). These tests demonstrate that our error 
model is correct, as the distributions very closely resemble a nor¬ 
mal distribution with zero-mean and unit variance shown by red 
curves on the panels. 


varez fc Plezll 19981 

et al.||2UU8[ ) model 


Plez||2Q12 ), the MARGS (Gustafsson 


vey line list (Rnffoni et al.|2Q14 Heiter et al.|2Q15 V5 for 


atmospheres and the Gaia-FbO bur- 

. , ^ ) T 

atoms and molecules). T'he grid includes ettective tem¬ 
peratures from 3000 K to 8000 K, and surface gravities 
from log^ = 0 to 5. Metallicities extend from as low 
as [Fe/H] = —5 with 1 dex steps until [Fe/H] = —3 and 
0.25 dex steps thereafter, extending past solar metallicity. 
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In the metallicity range applicable for this study, [a/Fe] 
ratios vary between 0.0 and +0.8 at steps of 0.2 dex. We 
redshift our interpolated spectrum by velocity V such 
that the normalised synthetic flux at an observed point A 
is given by S' (A [l + ^] , cj), where c is the speed of light. 
The observed continuum is modelled as a low-order poly¬ 
nomial with coefficients hj that enters multiplicatively: 


N-l . r 

M(A,aj,i;,{6}) = E ^channelj X S ( A 
3=0 ^ ^ 


i+’i 

C 


( 1 ) 

The continuum in each observed channel (HRIO and 
HR21) are modelled separately. In practice we found a 
first-order polynomial to sufficiently represent the contin¬ 
uum in each channel. Lastly, we convolve the model spec¬ 
trum with a Gaussian LSF (with free parameter TV) to 
match the resolving power in each channel, and resample 
the model spectrum to the observed pixels {A}. Although 
the spectral resolution TZ in each channel is reasonably 
well-known, recent refocusing of the GIRAFFE spectro¬ 
graph has improved the quoted spectral resolution. For 
this reason we chose to include the spectral resolution 
7^ as a nuisance parameter with reasonable priors and 
marginalise them away. The prior on spectral resolution 
was uniformly distributed to within ±30% of TZ =16200 
and 19800 for the HRIO and HR21 setups respectively. 
After convolution with the LSF, binning to the observed 
pixels {A} and assuming Gaussian error cr^, the proba¬ 
bility distribution p (F^|Ai, cr^, cj, V, {6}, {7^}) for the ob¬ 
served spectral flux Fi is: 


p{Fi\\i,(Ji,u,V,{b},{n}) 


1 


exp 


\ r 


( 2 ) 

Under the implied assumption that the data are inde¬ 
pendently drawn, the likelihood of observing the data 17, 
given our model, is found by the product of individual 
probabilities: 


N 


C = \{p{Fi\\uauUJ,V,{b),{JZ)) 


(3) 


and the probability 'P of observing the data is propor¬ 
tional up to a constant such that: 


V(xC{D\0) X Vr{0) 

\nV = \nC{D\e)F\nVr{e) (4) 

where Vr{0) is the prior probability on the model pa¬ 
rameters 0. 

In practice we found that the spectral range of our data 
were not particularly informative of the effective temper¬ 
ature Teff for very metal-poor stars. For these stars the 
data were prone to favour unphysically cool super giant 
stars of extremely low metallicity. We did not find the 
same effect for more metal-rich stars, where there are 
sufficient neutral and ionising transitions present to ac¬ 
curately constrain the stellar parameters. Given most of 
our candidates are indeed metal-poor, we found it pru¬ 
dent to fix the effective temperature using the DES pho¬ 


tometry and a colour-temperature relatiorj^ We found 
this had no significant impact on our posteriors for the 
foreground dwarfs - where the spectra are indeed in¬ 
formative of effective temperature - and ultimately did 
not substantially alter our inferred metallicity dispersion 
for either satellite. Only the mean satellite metallicities 
were affected. The uncertainties in effective temperature 
listed in Table were calculated by propagating the DES 
photometric uncertainties with the intrinsic uncertainty 
in the colour-temperature relation. Thus, our forward 
model is subject to our photometric temperatures and 
has only 10 parameters: log^, [Ee/H], [a/Ee], U, the 
resolving powers T^hrio and 7^hr21, as well as two con¬ 
tinuum coefficients in each channel. 

The initial model parameters V and uj (modulo Teff) 
were estimated by performing a coarse nor malisation of 
the data an d cross-correlating it against the |de Laverny| 
et al. (2012) grid. Although we have fixed Teff and previ¬ 
ously determined V (see Section 3.1) we still carried out 
the cross-correlation to yield a reliable initial estimate of 
log^, [Ee/H] and [a/Ee]. We also used the synthetic flux 
at the grid point with the peak cross-correlation coeffi¬ 
cient to subsequently estimate the normalisation coeffi¬ 
cients {b}. We numerically optimised the negative log- 
probability — In 7 ^ from the initial point using the Nelder- 
Mead algorithm (Nelder & Mead |1965). Eollowing opti¬ 
misation, we sanipled the resulting posterior using the 
affine-invariant Marko v Ghain Monte Garlo (MC MG) 


sampler introduced by Goodman fc Weare 


implemented by 4'oreman-lVlackey et al. 


C 

13b. In 


and 
all 

cases at least 200 walkers were used to explore the param¬ 
eter space for more than 2500 steps (> 5 x 10^ probability 
evaluations) to burn-in the sampler. These probability 
calls were discarded and the chains were reset before pro¬ 
duction sampling began. We tested our MGMG analyses 
for convergence by examining the auto-correlation times 
(e.g., ensuring high effective sample numbers per param¬ 
eter) and the mean acceptance fractions over time. We 
also re-ran a subset of our analyses with many more eval¬ 
uations (for both burn-in and production), verifying that 
there was no change to the resulting posteriors. 

We list the photometric effective temperatures and 
other inferred stellar parameters (given the model and 
effective temperature) in Table 1^ The middle panels of 
Eigure show the inferred sur^e gravities for all can¬ 
didates. Our confirmed members agree well with the 
metal-poor isochrone shown. This is particularly true 
for the higher-quality Reticulum 2 data, with the pos¬ 
sible exception of the metal-poor giant Reti 4, where 
the log^ seems quite low. Although this star has the 
highest S/N ratio in our sample, metal-poor super giant 
stars are very challenging to model from an astrophysi- 
cal perspective. Nevertheless, the marginalised posterior 
metallicity distribution for Reti 4 agrees excellently with 
stars further down the giant branch of lower S/N ratios. 
Table also lists the reduced values of our spectral 
fits. In general the values are quite close to 1. How¬ 
ever due to lower S/N and imperfect sky-subtraction in 
the Horologium 1 data, the y^/d.o.f. values for those 
candidates are slightly higher. 


^ The colour-temperature relation was defined by fitting the 
relation between DES g— r colors and effective temperature from 
SEGUE ( |Lee et al.|2008| . 
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Fig. 4. — Inferred radial velocities and stellar parameters for all Reticulum 2 (top) and Horologium 1 (bottom) candidates. Confirmed 
members (marked in red) cleanly separate from foreground contaminants in v — [Fe/H] space (left). Photometric effective temperat ures fsee 
text) and s pectroscopically-derived surface gravities are shown in the middle panel, overlaid upon a 12 Gyr PARSEC isochrone (|Bressan| 
|et al.|2012|) with a metallicity of [Fe/H] = —2. Histograms of the maximum a posteriori metallicities of confirmed members are shown m 
the right-Jiand panels. 


Our radial velocity determination for Reti 22 confirms 
it as a hori zontal branch member of Reticulum 2 (see 
Section 3.3). However, the photometric temperature es¬ 
timate of 8468^295 K prohibited us from inferring other 
stellar parameters for this star, as it is hotter than the 
boundary (8000 K) of the spectral grid. 

Two members of Reticulum 2 have maximum a poste¬ 
riori (MAP) log^ values that are consistent with being a 
dwarf (Figure®. This is inconsistent with the photom¬ 
etry, since the main-sequence is too faint for us to tar¬ 
get with standard exposure times for Gaia-ESO Survey 
Milky Way fields. However, the negative uncertainties on 
log g for these two Reticulum 2 members are considerably 
large, making them deviate from the giant branch by only 
1 to 1.5(7. As a test we constrained the prior on surface 
gravity to be uninformative between log^ G [—0.5,4.0], 
forcing the star to be a giant/sub-giant, but we found 
no statistically significant difference in the marginalized 
posterior metallicity distribution. 

Similarly while all confirmed Horologium 1 stars are gi¬ 
ants, and the foreground contaminants are clearly dwarfs, 
one star (Horo 17) has a very low S/N and consequently 
has an extremely large negative uncertainty in log^. 
While the posterior demonstrates the star is not a dwarf, 
we cannot place its precise location on the giant branch. 

Our inferred [a/Fe] abundance ratios are informative, 
even with their large uncertainties. Unsurprisingly, we 
found the [a/Fe] ratio to be strongly correlated with 
other stellar parameters, particularly [Fe/H]. We find 
the foreground contaminants to largely follow the well- 
studied Milky Way trend in [Fe/H] — [a/Fe]. All the 


TABLE 1 

Summary properties of Reticulum 2 and 
Horologium 1 



Reticulum 2 

Horologium 1 

a (J2000) [deg] 

53.9256 

43.8820 

d (J2000) [deg] 

-54.0592 

-54.1188 

Distance [kpc] 

30 

79 

Mv 

-2.7 ±0.1 

-3.4 ±0.1 

Ellipticity 

0 I5Q+0-02 

< 0.28 

ri /2 [arcmin] 

'^•'^^-0.12 

1 31 +0.19 
^•'^^-0.14 

’’1/2 [pc] 

32 IJ® 

30T^ 

Vhel [km s“l] 



a{V) [km s-i] 

q 00+1.64 
'^•^^-0.49 

^•^-0.9 

Mass(< ri/ 2 ) [Mq] 

2.35lo 13 10® 

5.5/:^+ X 10® 

M/Lv[Mq/Lq] 

479151 "^ 


[Fe/H] 

-2.46lo;i® 


f7([Fe/H]) [dex] 

'^•^^-0.05 

0 17+0-2 
' -0.03 

[a/Fe] 

0.40 ± 0.04 

0.30 ±0.07 


Note. — Prop erties above the horizo ntal separator 
were adopted from |Koposov et al.| (|2015|). Those below 
the separator were determined in this study. 


confirmed members in Reticulum 2 and Horologium 1 
appear to have at least [a/Fe] > +0.2, with an point esti¬ 
mate (assumed (i-function distribution) of [a/F^e] ^ +0.4 
in Reticulum 2 and ~ + 0.3 in Horologium 1. 

3.3. Joint modeling of satellites kinematies and 
ehemistry 

Having measured the chemical abundances and radial 
velocities of individual stars in two satellites the next 
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TABLE 2 

Positions, velocities, stellar parameters and membership for Reticulum 2 and Horologium 1 candidates 


Object 

a (J2000) 

(5 (J2000) 

9 

Vhel 

TeS 

logs 

[Fe/H] 

[a/Fe] 

Ared 

Member? 


[deg] 

[deg] 

[mag] 

[km s-1] 

[K] 





Reticulum 2 

Reti 0 

53.9424 

-54.1260 

19.27 

218.5 ± 1.1 

568012J^ 

q Qq + 0.35 
C).yc)_o 33 

— 1 99“^®'^^ 

J..»»_0 17 

0 55+®-!^ 

U.OO_o 19 

1.11 


Reti 1 

53.8133 

-54.1452 

19.78 

78.9 ± 1.8 

57141221 

C).yD_o 98 

9 q9+0.27 

^•'^^-0.25 

0 61+®'^^ 
’®-OJ-_0.22 

0.96 

Yes 

Reti 2 

53.8072 

-54.0824 

19.74 

60.0 ±2.1 

5729l2« 

-0.43 

—2 6‘i“*”®-^® 

Z.DO_o 32 

0.45l°;!! 

0.95 

Yes 

Reti 3 

53.9045 

-54.0670 

18.60 

65.6 ±0.9 

5558 I 19 I 


9 09+0.17 
^•°^-0.12 

0.65l°;}} 

1.00 

Yes 

Reti 4 

53.8494 

-54.0687 

16.47 

66.3 ±0.2 

4896i;“ 

u.Dy_o 02 

—2 50“*”®-®^ 
Z.OU_o 02 

0 23+®-®! 
O'-^'^-O.Ol 

1.43 

Yes 

Reti 5 

53.8374 

-54.0633 

18.97 

69.1 ± 1.0 

565512J6 

q q/yTO.63 
-0.26 

-2 54+®-18 

^•O^-0.16 

0 59+®'!^ 

U.oy_0 22 

1.03 

Yes 

Reti 6 

53.7260 

-54.0994 

18.97 

70.8 ± 1.1 

561712“ 

q ^0+0.69 

—2 56“*”®-^® 
Z.OO_o 25 

0 14+®-20 

0.91 

Yes 

Reti 7 

53.7399 

-54.0920 

18.97 

61.9 ±0.8 

55641201 

-0.17 

—2 03“*”®-^^ 

Z.U9_o 19 

0 39+®'!^ 

U.oy_0 15 

1.05 

Yes 

Reti 8 

53.7605 

-54.0650 

19.26 

65.4 ± 1.8 

566912“ 

q O^ + 0.89 

—2 51+®-^^ 
^•"^-^-0.26 

0 29+®‘^® 
'^•^^-0.22 

0.91 

Yes 

Reti 9 

53.8209 

-54.0675 

19.74 

62.9 ±3.7 

5671 

ODm_208 

4-25l°;“ 

3 1 8+®-30 
'J--L'J_o.45 

n 99+0.30 
'^•^^-0.20 

1.44 

Yes 

Reti 10 

53.8540 

-54.0418 

19.72 

65.6 ± 1.2 

5562“*”^®'^ 

OODZ_2 qq 

2.981°;“ 

1 88+®-^^ 
i.66_o 16 

0 42+®-^® 
'^•^^-0.22 

0.98 

Yes? 

Reti 11 

53.7986 

-54.0560 

19.35 

68.2 ± 1.7 

5651 

ODOi_2Q9 

^.»o_o 22 

—2 67“^®-^® 
-0.28 

0 39+®'^® 

u.oy_0 26 

1.17 

Yes 

Reti 12 

53.7896 

-54.0416 

18.27 

94.7 ±0.4 

567812“ 

5 07+0-01 
O-'O' -0.01 

n fTq+0.09 

O'- ' *^-0.09 

n 81+®-®’^ 

'O-'Jl-0.06 

1.02 


Reti 13 

54.0980 

-54.0885 

19.77 

220.3 ± 2.0 

56221201 

^•-'-O_o.90 

Z.UU_o 19 

0 58+®-!^ 

U.OC)_o 25 

1.04 


Reti 14 

54.0074 

-54.0681 

19.59 

63.4 ± 1.7 

565512“ 

^.»o_0 25 

—3 19+®-®! 

'^•-^^-0.33 

0 46+®-^® 
'O-^'J-0.27 

1.14 

Yes 

Reti 15 

53.9502 

-54.0638 

18.30 

63.5 ±0.5 

542ll“i 

9 71 +0.12 
^•'1-0.14 

— 1 98“^®'®® 
i.yo_0 05 

0 47+0-06 
-0.07 

1.11 

Yes 

Reti 16 

53.9582 

-54.0559 

19.71 

292.1 ± 1.3 

58281223 

4 84+0-15 
^•°^-0.28 

-l.39l°;“ 

0 54+®'!^ 
'0-O"±_o.23 

1.06 


Reti 17 

53.9845 

-54.0545 

18.90 

65.9 ± 1.2 

57241220 

3 49+0-26 

-2-68l°;^S 

0 22+®'^^ 
'O-^^-0.17 

1.09 

Yes 

Reti 18 

54.0323 

-54.0432 

17.46 

61.4 ±0.4 

5343l“3 

-0.06 

9 qr)+0.04 
^•'^^-0.14 

0 42+®-!® 
'^•^^-0.02 

1.05 

Yes 

Reti 19 

53.9923 

-54.0346 

19.32 

65.0 ± 1.4 

57411221 

0 91 +0.30 
'^•^-^-0.22 

9 q[r+0.23 

Z.C)0_Q 3Q 

0 22+®'^^ 
'^•^^-0.22 

1.06 

Yes 

Reti 20 

54.0163 

-54.0073 

17.79 

-23.6 ±0.3 

5 IO 7 IJ 7 I 

9:.DU_o 09 

-0-551^5 

0 19+®-®! 

1.06 


Reti 21 

54.0952 

-53.9987 

18.59 

128.9 ±0.5 

572712“ 


-0.88l°;°! 

—0 01+®-®® 
'O-’Ol-0.06 

0.97 


Reti 22 

54.0779 

-53.9625 

18.05 

61.6 ±2.6 

84681300 





Yes 

Reti 23 

53.8798 

-54.0300 

17.64 

59.6 ±0.5 

53861“! 

9 oq+0.09 
Z.OC)_o 09 

-2-681°;°! 

0 47+0-30 
-0.07 

1.16 

Yes 

Reti 24 

53.9127 

-53.9323 

17.88 

-19.2 ±0.3 

5138l“! 

4 48+®-®® 
^•^®-0.05 

-O-ll+ol 

0 09+®'®^ 

U.U^-O 02 

1.00 


Horologium 1 

Horo 0 

43.9692 

-54.3117 

18.81 

110.6 ±0.4 

51351 “0 

4.841°;}® 

-0-47l°;°! 

0 80+®'®® 
u.au-o 06 

2.01 


Horo 1 

44.0306 

-54.2768 

18.28 

254.2 ± 0.2 

4640l“! 

4-67l°;}2 

-0-88l°;“ 

0-29l°;°l 

1.58 


Horo 2 

43.8105 

-54.2267 

19.20 

152.6 ±0.8 

53651 “0 

4 95+0-05 

T.yO_o 12 

-0.59l°;}° 

0-38l°;}? 

2.23 


Horo 3 

43.6503 

-54.1802 

17.73 

27.5 ±0.1 

43891“! 

4 24+0.05 
^•^^-0.02 

—0 89“^®'®^ 
u.c)y_o 03 

0 21+®-®® 
'®-^i-o.oi 

2.03 


Horo 5 

44.1126 

-54.2174 

17.79 

64.8 ± 0.5 

50381}“ 

4 94+0-07 
^•^^-0.14 

— D 58+®-®® 

U.05_o 08 

0.29l°;°! 

1.58 


Horo 6 

44.1567 

-54.1941 

19.06 

91.2 ±0.9 

487ll“| 

4.79l°;}® 

—0 71+®-®® 

O'- '1-0.09 

0.38l°;“ 

2.22 


Horo 7 

44.0076 

-54.1986 

19.30 

163.1 ± 1.0 

51481 }32 

4.88l°;}° 

—0 46“^®'!^ 

0 89+®-!^ 
’^•'^^-0.10 

1.53 


Horo 9 

43.9179 

-54.1353 

18.71 

118.5 ±0.5 

49931}“ 

0 71+0.23 
O'-'1-0.15 

—2 55+®-ll 
Z.OO_o 18 

0 83+®-!^ 
^'•'^^-0.14 

1.97 

Yes 

Horo 10 

43.8967 

-54.1122 

19.31 

116.6 ±0.1 

45041}“ 

n 58+0-04 

U.OC)_o 02 

— 3 01“^®-®^ 
'^•'^1-0.03 

0 86+®'®® 
’®-'^®-0.03 

1.83 

Yes 

Horo 11 

43.9699 

-54.0877 

18.83 

114.6 ±0.7 

49721}“ 

1 80+0-30 
i.ou_0 29 

—2 79+®-!’^ 
'^-0.17 

0 83+®-!® 

U.40-0 18 

2.24 

Yes 

Horo 15 

43.8912 

-54.0939 

19.08 

105.6 ± 1.0 

50261}“ 

1 45+0-39 

l-40_o 41 

9 77+0.17 
' ' -0.22 

0.15l°;}} 

2.02 

Yes 

Horo 17 

43.8719 

-54.0727 

18.65 

108.1 ± 1.9 

52631 }3| 

8 1 [:;+l-63 
0--LO_0.55 

—2 86“^®'^^ 
Z.C)D_o 25 

0.18l°-+ 

2.03 

Yes 

Horo 18 

43.8497 

-54.0445 

18.40 

17.1 ±0.2 

49251}“ 

4 44^n'2n 

^•^^-0.09 

— D 58+®-®® 
U.05_o 05 

0 1 0+®-®5 

'®-l°-0.04 

4.25 


Horo 19 

43.8867 

-53.9968 

19.79 

18.7 ±0.6 

53031}“ 

A 48+0-27 
^•^®-0.20 

—0 84+®-!® 
O'-'^^-O.lO 

0.021°-“ 

1.57 


Horo 20 

43.6278 

-54.0217 

18.28 

98.5 ±0.2 

49541}“ 

4 03+®-^^ 
T.DC)_o 10 

—0 60+®-®® 
U.DU_o 06 

0-08l°;“ 

3.00 


Horo 21 

43.6917 

-53.9571 

19.36 

67.4 ± 1.1 

53421}“ 

4 95+®-!^ 

T.yo_o 17 

—0 12+®-^® 
'®-l^-0.14 

— 0 12+®'!^ 

2.96 


Horo 22 

43.9334 

-53.9473 

19.02 

70.0 ±0.6 

49791 }g! 

T.DC)_o 28 

—0 49+®-!® 

0 06+®'!^ 
U.UD-o 10 

1.99 


Horo 23 

43.8368 

-53.9240 

18.69 

44.2 ± 0.3 

48211}“ 

4 70+®-!^ 
^•'0-0.13 

—0 48+®-®’^ 

0 25+®-®® 
’®-^^-0.07 

1.73 
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step is combining all the available information in order 
to obtain the most reliable inference on the average ve¬ 
locity, metallicity and their dispersions, while properly 
accounting for any potential foreground contamination. 
Figure shows the radial velocity of stars versus the dis¬ 
tance to the center of two satellites. It clearly illustrates 
that although the velocity signal due to the satellites is 
quite prominent, the contamination - albeit minor - still 
has to be taken into account. Thus, in order to describe 
the velocity distribution of each satellite we adopt the 


following set of mixture m odels (see e.g. Walker et al. 
2009 Koposov et al.|[2MT for similar approach). 


P{V,^\4>)=fPsati^)J^{V\Vo,a) + 

il-f)P,gi^)mV\V,g,a,g) (S) 

Where V is the heliocentric velocity. A/’ is a Gaussian 
distribution, / is the fraction of objects belonging to the 
satellite, 0 is the shorthand notation for the parameters 
of the model and ^ are ancillary variables that help us 
identify members (such as metallicity and/or distance 
from the center of the object). The key assumption 
is that the RV distribution for each of the satellites is 
Gaussian and RVs of background/foreground stars are 
also Gaussian distributed (a reasonable assumption given 
very small number of such stars). 

Since each RV measurement comes with an error bar, 
the actual likelihood of each RV point V and error cji is 
the convolution of the model from Eq. with the Gaus¬ 
sian error: P{D\(j)) oc j P{V\(j)) J\f{V\Vi^cFi) dv. Given 
that the underlying velocity models P(V|0) are Gaus- 
sians themselves, the integral is trivial to analytically 
compute. 

The ancilliary parameters ip serve the purpose 
of helping to separate the satellite members from 
the background stars. For Reticulum 2 we use 
= [Fe/H] and model the joint distribution of 
metallicity and RV. We assume that the metallic- 
ities of both the background and the object are 
Gaussian-distributed (with different means and vari- 
ances): P.at([Fe/H]) = A/-([Fe/H]|[Fe/H]at), 
ng([Fe/H]) = A/'([Fe/H]|[Fe/H],g,a[pe/H].6g). For 

Horologium 1 we have fewer potential members, so we re¬ 
quire more information than RV and metallicity. There¬ 
fore we model the joint distribution of RV, metallic¬ 
ity, and distance from the center of the satellite: = 

{r, [Fe/H]}. The metallicities are modelled as Gaussian 
distributions, while an exponential density model with 
the morphological parameters from K15 is used to repre¬ 
sent the distance distribution of satellite member stars: 

= ( 6 ) 

where h is the exponential scale length. 

The model for the background sources assumes a uni¬ 
form distribution within the field Pbair) = 2rlr‘j, where 
the Vf is the field radius. 

The full list of parameters in our mem¬ 
bership modelling for both Reticulum 2 and 

Horologium 1 was /, Vat, ^sat, [Fe/H]^^^, Cr[Fe/H],sat 

and Vbg,abg, [Fe/H]^^, (T[Fe/H],6s, respectively. 

We adopt uninformative priors on Vsat^ Vhg^ 


and [Fe/H]^^, Jeffreys priors on the distribution disper¬ 
sions and /. The posterior was then sampled using the 
ensemble MGMG samp l er im plemented in Python by 


Foreman-Mackey et al.| (2013). The posteriors for the 
parameters of the satellites are shown on Figures and 
0 The resulting parameter measurements quoted in Ta¬ 
ble are the ID MAP values, and the uncertainties are 
the 68% percentiles. 

4. DISCUSSION AND CONCLUSIONS 
4.1. Reticulum 2 

On the basis of inferred kinematics and chemistry, our 
analysis has unambiguously identified 18 members in 
Reticulum 2. Of these stars, 17 are red giants and one is 
a horizontal branch star. 

We find an intrinsic velocity dispersion of 


3.22tl 


.64 
-0.49 


kms ^ in Reticulum 2. 


Although our 
value is slight l y lowe r tha n the parallel analy ses by 
Walker et al. (2015b) and Simon et al. (2015P the 
velocity dispersion measurements from all studies are 
consist ent within the uncer tainti es. As already pointed 
out by Simon et al. (2015) and Walker et al. (2015b), 
the velocity dispersion unambiguously indicate that 
Re ticulum 2 is a dwar f galaxy. Using the mass estimator 
of Wolf et al. (2010) we can estimate the total mass 
inside half-light radii to be 2.541q;32 x IO^Mq for 
Reticulum 2, which corresponds to a mass-to-light ratio 
of ^ 500. However the total mass and the mass-to-light 
ratio has to be treated with caution, as Reticulum 2 
is very elongated (axis ratio of 0.4) and is potentially 
being tidally disrupted (see K15), therefore the mass 
estimator could be significantly biased. 

We also find a substantial spread in overall r netallicity 
of cr([Fe/H]) = 0.29^0 05 ^®^* contrast to 


Simon et 


^ir](|2015), we have also inferred [a/Fe] abundance ratios 
for all satellite candidates observed through the Gaia- 
ESO Survey. Although the uncertainties on [a/Fe] are 
large for most of our confirmed members, we found the 
Reticulum 2 data tended towards high [a/Fe] ratios. In¬ 
deed, the lowest [a/Fe] ratio of our 18 confirmed mem¬ 
bers exceeds +0.2 dex. The resulting estimate of [a/Fe] 
for Reticulum 2 is [a/Fe] = 0.40 ± 0.04, consistent with 
observations C)f well-studied present-d ay dwarf galaxies 
( Tolstoy et al.|2QQ9| [Kirby et al.||2011 ). 

'I'here are slight discrepancies in the estimated mean 
metallicity of Reticu lu m 2 between t hi s stud y and 


Simon et al. (2015) 


Walker et al. (2015b) finds 
[Fe/H] = —2.67 Fq‘ 34, consistent with our measurement 

+0.09 
0 . 10 - 


Simon et al. 


of [Fe/H] = -2.46^ 
parable value of [Fe/H] = —2.65 + 0.07 fro m Ga H equiv- 


(2015) finds a com- 


alent widths. The quoted uncertainty by [Simon et al. 
(2015) are the lowest of all studies, but given our uncer¬ 
tainties, [Fe/H] = —2.65 is a mere 1.9cr deviation. We 
explored this possible discrepancy by searching the Gaia- 
ESO Survey for HRIO and HR21 spectra of HD 122563, a 
well-studied met al-poor giant star . HD 122563 is a Gaia 
benchmark star ([Jofre et al.[[2014[), and has atmospheric 
parameters comparable to the stars in Reticulum 2. The 
lowest S/N in any single exposure of HD 122563 was ^20. 


^ [Simon et al.[ ([20I5[) use the Gaia-ESO data as well in their 
analysis ot Keticiiliim 2 
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Fig. 5.— The measured RV and spatial distance from the satellite centroid for Reticulum 2 (left) and Horologium 1 (right). The bottom 
panels show the RV distributions. The dashed line on top panels indicates the half-light radius as measured by K15. Red circles are the 
stars with [Fe/H] < —1.5. In the case of Reticulum 2 the peak in the radial velocities at V~ 65 km s“^ due to the satellite stars is obvious. 
In the case of Horogolium 1 the peak at V~ 100 km s“^ is less prominent but still significant, because all the stars in that peak are within 
~ twice the half-light radius and have low metallicity as opposed to high metallicity background stars located at larger distances from the 
satellite center. 


We analy sed these data with the model described in Sec¬ 
tion |3.2[ except we found it necessary to use a fourth 
order polynomial to account for the continuum in the 
HD 122563 spectra. We find a MAP [Fe/H] = -2.79, 
in good agreem ent with the acce pted literature value of 
[Fe/H] = —2.64 (Jofre et al. 2014). If anything our metal¬ 
licity scale may be ^0.1 dex more metal-poor than the 
benchmark, t he op posite direction to the discrepancy in 


Simon et al. (2015). Nevertheless, our robust uncertain- 
ties make our measu reme nt reasonably cons istent with 


Walker et al.| (|2015b|) and|Simon et al. (|2015). 

An unexpected discovery was also rnade in the Gaia- 
ESO Survey Reticulum 2 data. The left panels of Fig¬ 
ure show the inferred radial velocity an d metallicity 
from the model described in Section 13.21 Two stars 
(Reti 0 and Reti 13) are present at Vhei ^ 219 kms“^ 
with indistinguishable metallicities of [Fe/H] ^ —2. The 
log^ is largely uninformative for these stars: the MAP 
values are consistent with a sub-giant star, but the 
uncertainties are sufficiently large that both a dwarf 
or sub-giant are equa lly p lausible. We check ed the 
Walker et al. (2015b) and Simon et al.| ( |2015 ) stud- 


ies for other s tars at compa^ r able ve locities! We found 
one match in Walker et al. (2015b) (star Ret2-153 in 
their nomenclature), which turned out t o be Reti 0, 
and u nsurprisingly both stars were in the [Simon et al. 


(2015) study, simply marked as ‘non-member s’ of Retic¬ 
ulum 2. A subsequent search in the surrounding Gaia- 
ESO Survey Milky Way field (e.g., non-Reticulum 2 
candidates that were in the same field) revealed a 
further two stars {a^5) = (53.69300,-54.17860) and 
(53.73709, —54.10720) with similar systemic velocities to 
Reti 0 and Reti 13 (223.7 kms“^, 221.3 kms“^). We 
have not inferred stellar parameters for these additional 
two stars. However, if we ignore the information that 
Reti 0 and Reti 13 have indistinguishable metallicities, 
a simple calculation using the average number of stars 
per kms“^ at RV ^ 200 kms“^ gives the significance 
of having 4 stars within ^5 kms“^ to be ^9 9.5% (af- 


Vitellsj 

2010 

). Analogous to the 300 km s“^ stream near 

Segue; 

hpa 

ha et al.|2009 Norris et al.|2010 Erebel et al. 


2013), this kinematic feature is yet another reminder of 
the h ighly substructured nature of the Milky Way ha lo 


(e.g., Schlaufman et al.||2009 Starkenburg et al.||2Q09 ). 


4.2. Horologium 1 

With only 5 confirmed members in Horologium 1, 
we are far more sensitive to stochastic sampling ef¬ 
fects than what we are for Reticulum 2. Nevertheless, 
we find a large kinematic dispersion of d.O^o.g km s“^ 

and a metallicity dispersion of 0.171 q o 3 , firmly grouping 
Horologium 1 with other known dwarf galaxies. Note 
that our posteriors on cr {V) and cr([Ee/H]) have con¬ 
siderable asymmetry towards higher velocity and metal¬ 
licity dispersions, which is primarily attributable to the 
low number of confirmed members in our sample. When 
more data becomes available, it is reasonable to ex¬ 
pect that a larger metallicity dispersion may be found 
for Horologium 1, as our MAP cr([Ee/H]) is the lowest 
reported measurement for comparable ultra-faint dwarf 
galaxies (Eigure [^. A lower metallicity dispersion is 
strongly disfavoured by our data, and would be incon¬ 
sistent with the large velocity dispersion we observe. We 
also find Horologium 1 to have [a/Ee] = 0.30 ± 0.07, 
consistent with the Milky Way dwarf galaxy population . 

According to the mass estimator of [Wolf et al. (2010) 
we estimate the total mass inside the half-light radii of 
Horologium 1 is 5.25lJ]^7g x 10 ^Mq (notice however a 
very big error bar). The mass-to-light ratio of ^600 is 
similar to that observed in Reticulum 2. 


4.3. Comparison to other ultra-faints 

Overall both Reticulum 2 and Horologium 1 systems 
seem to be quite representative samples of the other 
known ultra-faint dwarf galaxies. 

• The average metallicity of stars in both systems is 
very low (one of the lowest among dwarf galax¬ 
ies), but both dwarfs he well on the existing 
mass/luminosity - metallicity correlation (see bot¬ 
tom left panel of Eigure [^. 

• The metallicity spread, although uncertain, is sig¬ 
nificantly different from zero, matching what is ob- 
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Fig. 6 .— 2D and ID marginalized posteriors for the parameters of the chemo-dynamical modeling of Reticulum 2 


served in other dwarf galaxies (top left panel of 
Figure [^. It is possible though that the spread 
seen in Horologium 1 and Reticulum 2 is somewhat 
smaller than the spread of 0.5-0.7dex o bserved in 
other ult ra-faint systems such as Segue 1 (Simon et 
al.|2011| |, but this could simply be a result of small 
sample sizes. 


The dark matter content and the mass-to light ratio 
in the observed systems seem to agree well with the 
existing correlations with galaxy luminosity (right 
panel of Figure]^. This suggests that even with 
these new discoveries of ultra-faint dwarfs we have 
not reached the limiting density scale of dark mat¬ 
ter, which would infor m us about the elusi ve prop¬ 
erties of dark matter ( Gilmore et al.|[^QQ7 ). 


4.4. Possible association with the Magellanic Clouds 

Given their proximity to the LMG and the SMG on 
the sky, there exists an exhilarating possibility that some 
of the newly discovered satellites, including Reticulum 2 
and Horologium 1, have once been part of the Great Mag¬ 
ellanic Family. If such a connection proves true, there is 
hope to link the internal properties of the dwarfs (e.g., 
their dark matter content, the star-formation and the en¬ 
richment histories) with their orbital motion before and 
during the accretion onto the Milky Way. Thus, finally, 
an in-depth self-consistent picture of the UFS formation 
and evolution can be assembled. 

Gonsidering the total number of satellites discovered 
in the SDSS, VST ATLAS and PanSTARRs surveys, the 
relatively small patch of sky covered by the first year 
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Fig. 7.— 2D and ID marginalized posteriors for the parameters of the chemo-dynamical modeling of Horologium 1 


of DES observations appears unusually rich in satellites. 
According to K15, the over-density of satellites around 
the Magellanic Clouds is moderate but significant, with 
at least 3-4 objects possibly belonging to the LMC/SMC 
pair. Note however, that the above calculation does not 
account for the fact that for the faintest systems (e.g., di¬ 
rect analogs of Reticulum 2 and Horologium 1), the SDSS 
census is incomplete beyond 50 kpc. Therefore the num¬ 
ber of faintest dwarfs within the DES footprint - namely 
those with My > —4 - have to be estimated under the 
assumption of their Galacto-centric radial distribution. 
Given perfect freedom, it seems plausible to find a radial 
profile flat enough to produce as many faint satellites as 
have been discovered in the DES data. This, however, 
would seem to be in tension with the lack of discover¬ 
ies from surveys like VST ATLAS and PanSTARRs. In 


the absence of completeness estimates for the ongoing 
imaging surveys, we attempt to clarify the connection 
between the newly discovered satellites and the Magel¬ 
lanic Glouds by complementing their 3D positions with 
the radial velocity measurements obtained with the VLT. 

The satellites’ kinematics can be compared to predic¬ 
tions made from osmo logical zoom-in sim ulations. Eor 
example, according to Sales et al. (2011), the distribu¬ 
tion of the satellites in phase-space reveals the time of 
accretion of the Magellanic system. By finding one suit- 
able LMC analog in the high-resolution A quarius suite 
(Springel et al. 2008), Sales et al. (2011) convincingly 
demonstrate that a high concentration of the former 
LMC companions is expected in the Cloud’s vicinity if 
the LMC has only had one peri-center crossing. More 
recently, a systematic analysis of 25 LMC analogs in 
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Fig. 8. — Left panel. Mean metallic i tv (bottom panel) and me t allicit v spread (top panel) as a function of luminosity for faint stellar 
systems. The data from|McConnachie| (|2012|); |Willman &: Strader] (|2012|) were used. The measurements of Horologium 1 and Reticulum 2 
are shown by red circles. Uur measurements confirm that these recently discovered galaxies extend the existing mass metallicity trends, 
and show significant spreads in abundances similar to the ones observed among other dwarf galaxies and not in globular clusters. Right 
panel The dynamical mass inside half-light radius vs galaxy luminosity. Both Reticulum 2 and Horologium 1 are on the very faint end of 
the existing relationship. The dashed line shows the dynamical mass expected for the stellar population with M/L\/=2 (i.e., without dark 
matter). 


the ELVIS suite o f cosmolofflcal zoom -in simulations has 
been presented by Deason et al. (2015). Here, rather than 
a report on a case study, an evidence for a trend between 
the 2 : = 0 phase-space scatter of the LMC satellites’ and 


son et al. (2015), the observed distribution of satellites 
discovered m the DES data appears consistent with a re¬ 
cent, i.e. <2 Gyr, accretion. Eor such late events, the 
authors also provide a rough estimate of the total num¬ 
ber of past LMC satellites in the DES footprint: based 
on positions alone, there should b e at least 4 such ob jects 
in the current DES sample. Both Sales et al. (2011) and 


Deason et al. (2015) emphasize the role that kineniatics 
play in uncovering the origin of the Milky Way satellites: 
chance spatial alignments are possible, however these are 
in general less likely in the vicinity of the group’s central 
(e.g., in the LMC). 

The benefit of these N-body simulations is that they 
paint a fully consistent cosmological portrait of the Mag¬ 
ellanic Group, both in terms of accretion history as well 
as the amount of expected substructure. However, it 
is obvious that these simulations cannot match either 
the exact orbit of the LMC nor the presence of its mas¬ 
sive companion, the SMC. To complement the cosmo¬ 
logical N-body zoom-in runs, controlled simulations of 
LMC/SMC accretion can be mass-produc ed for a much 


larger range of the in-fall parameters (e.g., Nichols et al. 
|mi| . We will describe the outcome of such an experi- 
ment in the future (see Jethwa et ah, in prep.). Mean¬ 
while, we can shed some light onto possible links between 
the Magellanic Clouds and Reticulum 2 or Horologium 1 
by comparisons with the observed kinematics of the Mag¬ 
ellanic gaseous Stream (MS). The stream of neutral hy¬ 
drogen emanating from the Clouds has been mapped out 


the group in-fall time is presented. According to Dea- Besla et ^.[2010 


across tens of degrees, and is co mplemented with well- 
documented kinemat ics (see e.g., [Putman et al.||2003 
Nidever et al. 2008). Additionally, several nurnerica. 
models exist exp l aining the genesis o f the MS (see e.g. 


Diaz & Bekki 2012). 


Figure 1^ shows the positions of the two satellites in the 
space of MS longitude Tms, MS latitude Hms and radial 
velocity. The top left panel of the figure gives the loca¬ 
tions of Reticulum 2 and Horologium 1 with respect to 
the distribution of the column de nsity of H I gas in th e 
Magellanic stream as detected by|Nidever et al. (2008). 
We split the H I detections into three bins accoming to 
the latitude Hms, marked with red (0 ° < Hms < 14°), 
green (-9° < Hms < 0°) and blue (-23° < Hms < -9°) 
colors. The location of Reticulum 2 (Horologium 1) is 
also shown as blue (red) filled circle. The top right panel 
presents the false-RGB composite map of H I in the plane 
of the MS longitude Lms and h eliocentric radial velo city 
IAsr- As previously shown by Nidever et al. (2008), in 
this projection the Stream’s H I content torniis a broad 
band, typically ^100 kms“^ in extent. As indicated 
by the rapidly changing color, portions of the stream at 
varying Hms contribute different amounts of velocity sig¬ 
nal at given Lms- Note however that the steep velocity 
gradient as a function of the MS longitude is predom¬ 
inantly caused by the Solar reflex motion. This is con¬ 
firmed in the bottom left panel of Eigurej^ which displays 
the map of H I in the plane of Lms and Galacto-centric 
radial velocity Vqsr- The Vqsr signature of the stream 
remains largely constant up to Lms ^ —40° where it 
starts to decline, but slower than with respect to Vlsr- 
As obvious from both velocity maps described above, 
the MS H I is detected near the positions of Reticulum 2 
and Horologium 1 in the (Tms, V) plane. However, the 
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Fig. 9. — Compa rison of the kinematic s of Reticulum 2, Horologium 1 and the Magellanic Stream. Top Left: Column density of H I in the 
MS as detected by|Nidever et al.|(|2008|) projected in the MS coordinate system (Lms? ^ms)- For further exploration, the gas is split into 
three bins in latitude sJiowri m red (0° < Bms < 14°), green (—9° < Bms < 0°) and blue (—23° < Bms < —9°) color. The location 

of Reticulum 2 (Horologium 1) are shown by the blue (red) filled circle. Blue (red) dashed line marks the locations of the LMC halo chosen 
for the line-of-sight velocity predictions that are displayed in the bottom right panel. Top Right: Heliocentric radial velocity Rlsr as a 
function of the MS longitude Lms- The MS H I gas density in this phase-space projection is shown as a false-color RGB composite built 
with grey-scale density distributions from corresponding Bms bins shown in the left panel. Note the strong Rlsr velocity gradient, which is 
chiefiy caused by the Solar refiex motion. Both Reticulum 2 and Horologium 1 appear to lie close to the lower envelope of the MS velocity 
signal. Bottom Left: Same as the top right panel, but for Galactocentric radial velocity Vqsr as a function of the MS longitude Lms- Note 
that the transformation from the LSR to the GSR has significantly reduced the velocity gradient observed in the top right panel. While 
the velocities of Reticulum 2 and Horologium 1 are consistent with the lower range of the MS stream kinematics at the corresponding Lms 5 
the gas with comparable Vqsr is located mostly at Bms < 0°5 evidenced by the blue color of the region of the MS nearest to the satellites. 
This is further illustrated in the bottom right panel. Bottom Right: Motions of the MS gas with 0° < Bms < 14°. This map confirms that 
near Reticulum 2, the MS H I gas motions are different by approximately 100 kms“^. On the hand, near Horologium 1, there exists MS 
H I gas whose velocity is similar to that of the dwarf. The black solid line is the track of the LMG’s orbit (see K15). Dashed curves show 
the prediction for mean velocity of the LMG halo, i.e. the projection of the LMG’s velocity vector onto the line-of-sight. Black dashed 
curve corresponds to the lines of sight crossing the LMG halo at Lms = 0°? and the blue and red dashed curves correspond to the halo 
slices shown in the top and left panels. Reticulum 2 is ~ 80 kms“^ away from the halo prediction. However, Horologium 1 is a mere 
kms“^ away from the halo prediction. 


region of phase-space occupied by the dwarfs is domi¬ 
nated by the gas at negative Bms (he., around the SMC) 
as evidenced by the blue tint of this portion of the im¬ 
age. The bottom right panel of Figure illustrates the 
range of possible velocities of the MS with positive B 
only (i.e., those around Reticulum 2 and Horologium 1). 
At the longitude of Reticulum 2, there is a gap of or¬ 
der of 100 kms“^ between the grey-scale density map of 
the MS H I and the velocity of Reticulum 2. However 
near Horologium 1, the MS gas completely extends to 
the measured velocity of the satellite. 

The offset in the line-of-sight velocity between the trail¬ 
ing MS H I and the satellites is expected, as stripped gas 
can experience a drag force through interactions with the 
hot Galactic corona. The drag will decelerate the stream 
clouds, causing them to fall to lower Galactocentric radii. 
To make sense of the radial velocities of Reticulum 2 and 
Horologium 1, let us recall their locations with respect 
to the Milky Way center and the Magellanic Clouds. In 
3D space. Reticulum 2 is in front of the Clouds, ^24 kpc 
away from the LMC and ^39 kpc from the SMC, while 
Horologium 1 is behind the Clouds, ^38 kpc away from 
the LMC and ^32 kpc from the SMC. Importantly, both 
satellites are trailing the LMC, as evidenced from the 


Cloud’s orbital motion shown in Figure 20 of K15. 

The black solid line in the bottom right panel of Fig- 
urej^shows the projection of a backward-integrated LMC 
orbit from K15, namely the one with the NFW’s concen¬ 
tration c = 10, projected onto the plane of Lms and 
kcsR- The orbit attains negative line-of-sight velocities 
at Lms < —20 ° and appears to be in reasonable agree¬ 
ment with the velocity of Horologium 1 at the corre¬ 
sponding MS longitude. This implies that the 4D coor¬ 
dinate of Horologium 1 is consistent with those expected 
for the LMC’s trailing debris. At the location of Retic¬ 
ulum 2 the velocity gap of ^100 kms“^ persists. It is 
tempting to assert that this precludes the possibility of 
an association between Reticulum 2 and the LMC. How¬ 
ever, it is unrealistic to expect Reticulum 2 to behave 
simply like the LMC trailing debris. During tidal dis¬ 
ruption it is normal to expect the trailing debris to form 
from particles with higher energy and angular momen¬ 
tum than the progenitor. Subsequently the trailing de¬ 
bris also have larger Galactocentric radii on average, and 
longer orbital periods. Yet, we know that Reticulum 2 
is closer to the Milky Way’s center than the LMC itself. 
Thus, in order to explain its origin as part of the Mag¬ 
ellanic family, an additional factor needs to be included. 
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It is conceivable that an interaction with the SMC would 
be sufficient to drive Reticulum 2 onto it’s current orbit. 

Finally, we have not yet considered the possibility that 
Reticulum 2 and/or Horologium 1 could still be bound 
to the LMC. Superficially, such situation seems unlikely 
given the dista nces between the LM C and the two satel¬ 
lites. However, Munoz et al.|(2006) report spectroscopi¬ 
cally confirmed detection ot the likely LMC’s stellar halo 
some 22° degrees away from the LMC, at angular dis¬ 
tances comparable to Reticulum 2 and Horologium 1. 
Motivated by this discovery, we test whether the line- 
of-sight velocities of Reticulum 2 and Horologium 1 are 
consistent with the LMC’s halo. The dashed lines in 
the bottom right panel of Figure show projections of 
the LMC velocity vector onto several lines-of-sight. The 
black dashed curve corresponds to the lines-of-sight slic¬ 
ing the LMC’s halo at Lms = Blue (red) dashed curve 
shows the run of the projection of the mean LMC’s halo 
velocity along the line-of-sight moving from the LMC’s 
center to the location of Reticulum 2 (Horologium 1) as 
shown in the top left panel of the figure. At the po¬ 
sition of Reticulum 2, the line-of-sight velocity of non¬ 
rotating LMC’s halo would be ^ — 20 kms“^, some 
^80 kms“^ away the Galactocentric velocity of Retic¬ 
ulum 2, ^ —100 kms“^. However, at the position 

of Horologium 1, the LMC’s halo’s line-of-sight veloc¬ 
ity is predicted to be ^ —50 kms“^, only 15 kms“^ 
away from the Galactocentric velocity of the satellite 
^ kms“^. It is important to note that while 
Horologium 1 is located at a similar angular distance 
from the LM G as compared to the stellar halo detections 
reported by Mufioz et al. (2006), it is probably twice 
as far away from the LlVlG in 3D, i.e. 40kpc instead 
of 20kpc. To explain the dynamics of the LMC+SMC 
system and the formation of the MS, |Besla et al. (2010) 
advocate the existence of a massive (Mvir = 3 x 
dark matter halo of the LMG. The corresponding virial 


radius of this halo would be > 100 kpc, while its (ap¬ 
proximate) tidal radius just under 40 kpc. If the LMG is 
on its first pericenter crossing, then given the weak tides 
the LMG has been experiencing so far, it is reasonable 
to expect that its dark matter distribution can extend 
as far as 40 kpc. Therefore, surprisingly, there appears 
to be some probability that Horologium 1 is still gravi¬ 
tationally bound to the Large Magellanic Cloud. 


Based on data products from observations made with 
ESO Telescopes at the La Silla Paranal Observatory un¬ 
der programme ID 188.B-3002. These data products 
have been processed by the Cambridge Astronomy Sur¬ 
vey Unit (CASU) at the Institute of Astronomy, Univer¬ 
sity of Cambridge, and by the FLAMES/UVES reduc¬ 
tion team at INAE/Osservatorio Astrofisico di Arcetri. 
These data have been obtained from the Gaia-ESO Sur¬ 
vey Data Archive, prepared and hosted by the Wide Eield 
Astronomy Unit, Institute for Astronomy, University of 
Edinburgh, which is funded by the UK Science and Tech¬ 
nology Eacilities Gouncil. 
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